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ERROR ANALYSIS FOR RELAY TYPE 
SATELLITE-AIDED SEARCH AND RESCUE SYSTEMS 

John W. Marini 


ABSTRACT 


An analysis is made of the errors in the determination of the 
position of an emergency transmitter in a satellite-aided search 
and rescue system. The satellite is assumed to be at a height 
of 820 km in a near circular near polar orbit. Short data spans 
of f'’>ur minutes or less are used. The error sources considered 
arc measurement noise, transmitter frequency drift, ionospheric 
effects and error in the assumed height of the transmitter. The 
errors are calculated for several different transmitter positions, 
data rates and data spans. The only transmitter frequency used 
was 406 MHz, but the results can be scaled to different frequencies. 

In a typical case, in which four Doppler measurements were taken 
over a span of two minutes, the position error was about 1.2 km. 
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ERROR ANALYSIS FOR RELAY TYPE 
SATELLITE- AIDED SEARCH AND RESCUE SYSTEMS 

INTRODUCTION 

In the operation of a satellite-aided search and rescue system, distress signals 
from an emergency transmitter are received by an orbiting satellite. The posi- 
tion of the transmitter is determined by an analysis of the Doppler shift of the 
fi - ^quency of the received signal. 

A practical consideration is the way by which the Doppler information is trans- 
mitted back to the ground for use by a rescue team. One method is to measure 
and record the Doppler-shifted signal at the satellite. The record is later 
"dumped" when the satellite passes within view of a system ground station [ 1) . 
Another method involves no recording but instead requires the existence of a 
number of ground stations so positioned that the satellite has at least one of 
them in view when it receives the emergency transmission, hi this case the 
Doppler Information can be immediately relayed to ground station. This can be 
done either by measuring the Doppler shift at the satellite and encoding the 
measurement onto the transmission from the satellite or by merely routing the 
emergency signal through a transponder [1]. Both the recording and the 
immediate relay method may be employed in the joint U. S.-Canada tests using 
the TIROS-N satellite (2, 3J. 

In a search and rescue system employing the relay method, Doppler data is 
obtained only during the interval when the satellite is in view not only of the 
emergency transmitter but also of the receiving ground station. Since the 
number of ground stations is limited by considerations of cost, the duta- 
collectlon Interval may be limited to only a small portion of the satellite pass 
over the emergency transmitter. It therefor becomes important to determine 
how the accuracy of the Doppler position determination is affected by the use 
of short data spans. The case where data is available from the entire satellite 
pass has been treated by Koch [ 4 ]. 

A second practical consideration is the rate at which Doppler information is 
transmitted. Both because the satellite might be requii'ed to handle a large 
number of emergency signals simultaneously and also in order to reduce the 
average power requirements of the emergency transmitters, it is advantageous 
to have the tiansmitters operate not continuously but instead repetitively in 
short bursts. The selection of the repetition rate involves a trade-off, since 
the improvement in the signal capacity of the 'satellite and the reduction of the 
power requirement of the emergency transmit er obtained by the u:<e of a low 
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repetition rate must be weighed against the error introduced in the calculation 
of transmitter positions by the reduction of the number of data points. 

In this report the error in the position calculation will be estimated using data 
spans of several durations and repetition rates for a number of transmitter 
locations with respect to the satellite ground track. The error sources 
considered will be Doppler measurement noise, assumed to be independently 
distributed for each point, linear drift of the emergency transmitter frequency, 
ionospheric effects as modeled in Reference 1, and error in the assumed height 
of the emergency transmitter above sea level. The latter error exists because 
it is usiuily better to estimate the transmitter height rather than solve for it 
when only a small quantity of Doppler data is available. 

The effect of satellite ephemeris errors is not treated. It is obvious, however, 
that both along-track and cross-track satellite position errors that remain 
nearly constant over the data span will result in an approximate one-to-one 
displacement of the calculated transmitter position. A constant radial (height) 
error in satellite position would have roughly the same effect as a corresponding 
negative error in the emergency transmitter height assumed. 

The ambiguity problem [5], the existence of an extraneous solution to the least 
square equations for transmitter location, is also not treated. In the mathematical 
method employed, position location errors are calculated under the assumption 
that the proper least squares solution has been determined. 

MATHEMATICAL METHOD 


The Doppler shift of the received signal is the difference between its frequency 
f, and that of the transmitted signal fj. Ignoring atmospheric effects and re- 
ceiver noise, the classical relationship between the Doppler shift of an uplink 
signal and the range rate p is 

P = -C(f,-f,)/f, (1) 

where c is the speed of light. In the satellite-aided search and rescue system 
under consideration, however, the exact frequency of the emergency transmitter 
ordinarily will not be known. Consequently a nominal frequency f„ must be used 
in (1) to calculate the range rate. Defining the range-rate so calculated to be 
the observed range rate 
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Taking the difference between (1) and (2) 


P-Po = -c 


Vn 


(3) 


Since the fractional variation in the received signal is small, the right hand 
side of (3) may be approximated* by a constant provided that the nominal fre- 
quency used is close to the transmitted frequency and that the latter is suf- 
ficiently stable over the duration of the satellite pass. Consequently the range 
rate observed at the time of the i'th Doppler measurement is modeled as the 
sum 


Poi = Pj + b (4) 

where p. is the true range rate and b is an unknown fixed bias arising from the 
use of a nominal transmitter frequency in the calculation of p^.. The relation- 
ships (4) hold only approximately because of measurement noise, neglect of 
atmospheric effects on propagation, and other modeling errors. 

The procedure usually employed to solve the family of equations (4) for the 
unknown transmitter position is the method of least squares. Since the satellite 
ephemeris is assumed known, the range rate Pj between the satellite and a given 
transmitter position becomes a function of transmitter position alone 


Pj = P^(<P, X, h) (5) 

where 0 and X are the transmitter latitude and longitude respectively, and h is the 
hei^t above the surface of an ellipsoidal Earth. In a single pass solution, the 
height of the transmitter is usually not solved for but histead is taken to be some 
assumed value h,. The least squares solution of (4) consists, then, of the values 
0 , X, and b that minimize the sum of squares 

S = S [p^,-p.(0,x,h^)-bl2 (6) 

i 

*This approximation, which ii often used [ J. 4) simplifies the least squares equations for determining the 
position of the transmitter. It is not a requirement, however, since the more accurate equation (1 ■ p^) f^ ■ f^, where f^ 
is treated as an unknown constant, can also be solved by non-Unear least squares. 
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The method used here to determine the magnitudes of the errors 

A0 = 0 ~ 0 

AX = X - X (7) 

Ab = b — b 

does not employ actual observed data, of course. Instead the position coordinates 
0, X , h are chosen, and the observed rango rate is taken to be the sum of the 
true range rate p (0, X« h), an error term e. , and optionally a bias term b. 


Poi = ^i(^’^»h) + b + €j (8) 

Substituting (8) into (6)» and linearly ejqianding p^ (0, X > h, ) about the values 
0) Xf h the sum of squares (6) becomes approximately 

r 3pj dpj 3p y 

The minimization of (9) is a standard problem in linear least squares. Using 
the matrix notation 


30 3X 

3pj 3pj 

30 dX 


30 3X 



( 10 ) 


( 11 ) 
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+ (h-hj) 

3p,/dh 

dpj/ 3h 

• 

• 

• 


* 

• 

3p„/ 3h 


(12) 


where n is the number of Doppler measurements, the minimizing values of A0, 
AX, and Ab are given by 


X = (A‘*‘A)-' a'^y, 


(13) 


being the transpose of the matrix A. 

Equation (13) is used to calculate the errors caused by oscillator drift, neglect of 
ionospheric effects, and error in assumed height. To determine the effect of a 
constant transmitted frequency drift df , /dt , for example, the components of Y 
are set equal to 


Y 


i 


(df,/dt) 


(14) 


and equation (14) is solved for X for each selected tramsmitter position (<*>, X* h). 
Similarly ionospheric effects are calculated by setting 


Y, = AR, 


(15) 


where AR| is the ionospheric correction given in Reference 1. An electron 
content of 6.5 X 10*^ electrons per square meter was used. 

The effect of a 1 kilometer error in the assumed height h, is calculated from the 
value of Y obtal'^ed by setting h - h, « i km and Cj » 0 in (12). 
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To determine the effect of measurement noise on and one takes the outer 
matrix product of X with Itself, From equation (13) 


XX’’ = (A'^A)*‘ A'*'YY'*‘A(A'^A)-* (16) 


The range-rate error is related to frequency error by the equation 


«i = -<c/f„)£n 


(17) 


where e^, is the error in the i'th measurement of the Doppler frequency. The 
errors are assumed to be independent with mean zero and a standard deviation 
Of . It follows that the range-rate noise errors 6^ are independently distributed 
with zero mean and the standard deviation a 


a = (c/f„)Of 


(18) 


Hence the expected value of the product YY^ that appears in (16) is simply the 
identity matrix times the variance o^. Taking the expected value of both sides 
of (16), therefore, gives 


EIXX^] = (AA)-‘o2 


(19) 


Consequently the covariance matrix of the position and bias errors per unit 
variance in range rate is (AA)~^ , and the standard deviations of A^and AX per 
unit standard deviation in range-rate noise are given respectively by the square 
roots of the first and second diagonal elements of the matrix (A A)*‘. 

SIMULATION GEOMETRY 

A circular orbit, similar to that planned for TIROS-N, with a radius of 7200 
kilometers and an inclination of 96. 7* was used in the simulations. The exact 
latitudes of the transmitter positions had little effect on the results. All of 
the positions used had latitudes of about 34*N. 
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Three transmitter positions were used to obtain the results given in Tables 1-4. 
These positions were taken so thr.t the satellite elevations as seen by an observer 
located at the transmitters were 20®, 50", and 80" at the time of closest approach 
(TCA) of the satellite to the transmitter. 

RESULTS 

The Doppler data used in the error analysis consisted of points taken at equal 
intervals over intervals of 1, 2 , or 4 minutes. In Table 1, an example is given 
at a 2 minute data interval consisting of seven Doppler points. The interval is 
centered at two minutes after TCA, Tables 2-4 show the results of the error 
analysis for each of the three stations. In Figures 1-4 the error caused by 
measurement noise is plotted as a function of elevation angle, number of points, 
length of the data interval, and data center displacement respectively. In the 
case of a large number of points, the noise error should be inversely proportional 
to the square root of that number, as shown in Figure 2. 

In Figure 5, height error is shown as a function of elevation angle. The curve 
is rou^ly approximated by the tangent of the elevation angle. 

The transmitter frequency was 406 MHz in all cases, but the results are easily 
scaled to other frequencies. Height error is independent of frequency. Noise 
and drift errors are scaled to a new frequency f by multiplying by the factor 
(406 /f), as follows from equations (17) and (14). Ionospheric effects vary as 
the inverse square of the frequency, the scaling factor being (406/f)^. They 
also vary directly with the integrated electron density used. Ionospheric error 
can be reduced, of course, by incorporating an ionospheric model in the least 
squares solution. 

The largest number of points used in a data interval was seven. Drift, height, 
and ionospheric errors show little variation with the number of data points per 
interval. Noise error for a larger number of points n may be estimated by 
multiplying the result for seven points by the factor V 7/n. 

As a numerical example, consider a 406 MHz transmitter whose drift rate is 
about 10 Hertz over 10 minutes. Assume that the transmitter is known to be at 
sea level, and that the standard deviation of the received signal is 0.29 Hz, 
corresponding to a determination that the frequency falls with uniform probability 
within a one Hertz band, the maximum elevation angle is 20"; and if four 
points are measured over a two minute interval centered at the time of closest 
approach, then the noise, drift, and height errors are 0.78, 0.87, and 0 
respectively. The ionospheric error varies directly with the electron content, 
and would probably be less than 0.32 km, which can be Ignored. The expected 
value of the error is estimated by taking the square root of the sum of the 
squares of the individual errors, giving about 1.2 km. 
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An analysis has been made of the errors in the determination of the position of an 
emergency transmitter in a satellite-aided search and rescue system. The 
analysis assumed a satellite at a height of 820 km in a near circular near polar 
orbit. Short data intervals of four minutes or less were used. Calculations were 
carried out assuming a 406 MHz signal, but the results can easily be scaled to 
other frequencies. Enough cases were treated to make possible at least a rough 
interpolation of the results to different data rates and different transmitter 
locations. 
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Table 1 

Example of Data Interval, Length: 4 min, Center: TCA + 2, 7 Points, 
Satellite Elevation: 20° at TCA. 


Time 

Elevation 

Range Rate 

9» 40” 

20° 

0 Km/sec 

10* 20” 

19.6° 

1. 05 Km/sec 

11* 

18.3° 

2. 02 Km/sec 

11' 40” 

16.5 

2. 88 Km/sec 

12* 20*' 

14.3° 

3. 59 Km /sec 

13* 

11. 9° 

4. 17 Km/sec 

13* 40** 

9. 5° 

4, 63 Km/sec 
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Table 3 

50® Satellite Elevation At Time of Closest Approach (TCA) 




m-* \Q \0 ^ 

CN ' OO 00 to i/^ tJ- o On On O *0 ro no 

O0OOOO(N<NfS0OOOOO<N<NfN0 00 0 


ONONO^^^oooooQC'ir-irNioootofncN'^ 

rnrOrfrorOfOroror^TtTt^mrorof^rorOTfTf'^ro 


ooor-r^<NO 

SSS8S^^^?i?SSSSSg'-$55§S8S 
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00 fN On Tt to t^Tfr^ON 

—^ONt^rOfNt^roOvDvOsDtofNOO^t^fNr^f^rocjoo 

r^-^-^fs-^ONod-^^dddfOfNoioNr^^ocdd 


4 0.72 0.050 1.35 0.26 

7 0.61 0.046 1.35 0.26 

TCA + 4 3 33.4 6.2 5.0 0.41 1.30 2.1 

4 3.9 0.35 1.30 1.70 

7 3.1 0.33 1.30 1.53 
















80“ Satellite Elevation At Time of Closest Approach (TCA) 

Station Latitude = 34, 2° Frequency g 406 MHz 

DATA INTERVAL | SATELLITE EL (deg) j RSS POSITION LOCATION ERROR (Km) 







TCA + 4 3 39.9 { 6.8 19.3 1.86 

4 ! 14.1 1.55 

7 i 10.8 1.40 
























DATA INTERVAL: 

1 MINUTE, THREE POINTS, CENTERED AT TCA 






NOISE ERROR (Km/Hz) 






ELEVATION - 50° 

CENTER OF DATA INTERVAL > TCA 

CURVE A NO. DATA POINTS ° 3 

CURVE B CONSTANT REPETITION RATE, 

1 POINT EACH HALF MINUTE 



1.5 2 2.6 3 4 6 6 7 8 9 10 

LENGTH OF DATA INTERVAL (minutes) 

Figure 3. Noise Error vs. Length of Data Interval 
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ELEVATION - 50“ 

DATA INTERVAL ^ 1 MINUTE 
THREE POINTS 



2 

DATA CENTER DISPLACEMENT FROM TCA (minutes) 

Figure 4. Noise Error vs. Data Center Displacement 
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